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Abstract 

We present work in progress on employing domain wall fermions to 
simulate N = 1 supersymmetric Yang-Mills theories on the lattice in 
d = 4 and d = 3 dimensions. The geometrical nature of domain wall 
fermions gives simple insights into how to construct these theories. 
We also discuss the obstacles associated with simulating the N = 2 
theory in d = 4. 



1 Chirality and accidental supersymmetry 

There has been intense interest in supersymmetry (SUSY) in the past two 
decades. The past several years have witnessed many interesting and com- 
pelling speculations about strongly coupled SUSY theories. It would be 
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interesting to test these conjectures on the lattice. However, since Poincare 
symmetry does not exist on the lattice, supersymmetry does not either. In 
principle, one could tune lattice theories to the SUSY critical point. How- 
ever, just as Poincare symmetry is recovered without fine tuning, one might 
hope that SUSY could be similarly obtained in the continuum limit. 

The secret to why the Poincare symmetric point takes no work to find 
is that the imposition of hypercubic symmetry and gauge symmetry ensures 
that Poincare symmetry is an accidental symmetry. All allowed operators 
that violate the symmetry are irrelevant. Therefore the continuum limit of 
the theory automatically exhibits more symmetry than it possesses at finite 
lattice spacing. If supersymmetry could arise as an accidental symmetry as 
well, then simulation of such theories would not entail fine-tuning of param- 
eters, and would be relatively simple. 

The outlook for this approach in SUSY theories with scalar fields is 
poor. ..scalar mass terms violate SUSY, are relevant, and cannot be forbidden 
by any symmetry, unless the scalars are Goldstone bosons. I will return to 
this issue later, when we talk about N = 2 super Yang-Mills (SYM) theories. 

However, there are some SUSY theories which do not entail scalars. Of 
particular interest is N = 1 SYM theory in d = 4 dimensions. In this theory, 
the only relevant SUSY violating parameter is the gaugino mass, which can 
be forbidden by a discrete chiral symmetry. Thus if one can realize the chiral 
symmetry on the lattice, N = 1 SUSY can arise as an accidental symmetry 
in the continuum limit. This is where domain wall fermions [JTJ, [| come 
in, for which chiral symmetry violation (for weak coupling) tends to zero 
exponentially fast in the domain wall separation. In this talk we clarify how 
domain wall fermions may be used to study SYM theories []. Throughout this 
talk we will actually discuss only the continuum version, as it is simpler to 
formulate, if less rigorous, and there are no technical or conceptual obstacles 
to translating this work to the lattice. We address in turn N = 1 in d = 4, 
N = 1 in d = 3, and N = 2 in d = 4. 



1 It has long been recognized that SYM theory can arise accidentally as the low energy 
limit of a theory with gauge and chiral symmetry, and the correct fermion representation 
||. Lattice implementation of iV = 1 SYM with Wilson fermions and fine-tuning is 
discussed in [1,0, [| [t| . Using domain wall fermions for simulation of N = 1 SYM was 
suggested in |£ , |)[ 10), but here we follow a different approach. Our results parallel prior 
work on N = 1 SYM theories in the overlap formulation jy], [f2|], which is equivalent to 
domain walls with infinite separation. 
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2 N = 1 SUSY Yang-Mills theory in d = 4 
dimensions 



N — 1 SUSY Yang-Mills theory in d = 4 Minkowski space consists of a 
gauge group with a massless adjoint Majorana fermion, the gaugino. It is 
expected to exhibit all sorts of fascinating features, such as confinement, 
discrete vacua, domain walls, and excitations on these domain walls which 
transform as fundamentals under the gauge group [p~3|] . 

The gaugino should arise as an edge state in a 5-d theory of domain wall 
fermions. The only subtlety in using the machinery of domain wall fermions 
is how to obtain a single Majorana fermion in Minkowski space, since without 
modification, the theory gives rise to massless Dirac fermions in Euclidian 
space. 

Let us first review how a massless Dirac fermion arises in the domain wall 
approach. Consider a Dirac fermion in a 5-dimensional Euclidian continuum, 
where the fifth dimension is compact: 25 = R8, 9 G (— ir, ir]. The mass of 
the fermion is given by a periodic step function 

m ( X5 ) = Me(6) = [_ M oth / erwise " (1) 

We introduce gauge fields independent of the coordinate x 5 , so that the 
Euclidian action is given by 

S 5 = J d 5 x i^D(x 5 )y , D(x 5 ) = [p 4 + d 5l5 + m(x 5 )] , Tm = T J . (2) 
Here If) ^ is the usual d = 4 gauge covariant derivative for a Dirac fermion in 



the adjoint representation. It is convenient to expand Psi and ^ as [14 



M.Xft, X 5 ) = En [b n (x 5 )P + + f n (x S )P-] ^n(<fy) , 
^(X^,X 5 ) = EnV'n(^) [b n (x 5 )P- +f n (x B )P + ) . 

Here P± — (1 ± 7s)/2 are the chiral projection operators, ip n and i\} n are 
ordinary 4-d Dirac spinors, and b n , f n form a complete basis of periodic 
functions satisfying the eigenvalue equations 

[d 5 + m(x 5 )}b n = fafn , [~d 5 + m(x 5 )]f n = [i n b n . (3) 

With this expansion, the action S5 may be rewritten as a theory of an infinite 
number of 4-d flavors with masses /i n , 

S 5 = J2j d ^ X [* ^4 + *A*n] i>n(x) ■ (4) 
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Figure 1: The left- and right-handed zeromode components, localized at the 
mass kinks. 

It is straight forward to solve the above equations for fi n . First of all, one 
finds zero modes 

/io = 0, b (x 5 )=e~r m ^ d y , f Q ( x ,) = e + P m ^ dy . (5) 

Note that bo is localized at 9 = — 7r/2, while /o is localized at 9 = +ir/2. 
Nonzero modes have wave functions which are linear combinations of sines 
and cosines appropriately matched at the locations of the domain walls. The 
corresponding eigenvalues are doubly degenerate 

fl n = y/M 2 + n 2 /R 2 , n= ±1, ±2,... (6) 

If instead of having a kink-like mass profile for the \& fermions we had a 
constant mass M (again with periodic boundary conditions), the correspond- 
ing eigenvalues ft would be 

^ = M, fi n = ^/M 2 + n 2 /R 2 , n = ±l,±2,... (7) 

Note that for n ^ 0, the eigenvalues //„ and equal. It follows that 

the ratio of fermion determinants for a kink and a constant mass is given by 
(assuming appropriate regularization) 

det[i(p A + 7 5 d 5 + Me(8))] det [ip 4 ] 

det [z( p 4 + 75 9 5 + M)\ det [i{ p A + M)\ U 
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Note that the right hand side of the above equation corresponds to a massless 
Dirac fermion and an uninteresting Pauli-Villars field. The left and right 
handed components of the massless Dirac fermion correspond to the edge 
states 6q and fo (see Fig. 1.). This method for obtaining a single massless 



Dirac fermion is robust when transcribed on the lattice [15|| : the beauty of the 
method is that there is no chirality in 5-d, and if one shifts or renormalizes 
the fermion mass term in the 5-d theory by 5m (with \Sm\ < M), the effective 
4-d theory still has a massless mode. That is because there is a gap in the 
bulk, so that 60 and fo fall off exponentially, while a chiral symmetry breaking 
fermion mass must be proportional to the (exponentially small) overlap of 60 
and Jo- 
in order to simulate N = 1 SYM theory, we need to impose a Majo- 
rana condition on ipQ. Note that in 4-d Minkowski space, the Majorana 
condition is ip = Cip T , where C is the charge conjugation matrix, satisfying 
(7 _1 7 M C = — 7^ T and C~ x T a C = —Tj for generators T a of real or pseudo-real 
representations of the gauge group. In Minkowski space charge conjugation 
interchanges left- and right-handed particles. In our Euclidian domain wall 
theory, the left- and right-handed modes live on the two different kinks. This 
suggests that the correct "Majorana" condition for the 5-d Euclidian theory 
is to define a 5-d reflection which interchanges the two chiral zeromodes, 
7Z 5 : 9 — > —6, and to impose the constraint on the 5-d Dirac fermions 

* = TZ 5 CW (9) 

The 5-d path integral then results in a fermion pfaffian, rather than a fermion 
determinant: 



Z5 = Pf[zft 5 C(f 4 + 75<95 + m0r 5 ))] • (10) 

It is straightforward to check that we are taking the pfaffian of an antisym- 
metric operator, as is required []. In terms of the mode expansion in 4-d 
fields, note that 1Z$ interchanges b n {x^) <-> /n^s), so the constraint yields 

= En[bn(x 5 )P + + f n (x 5 )Pj\ C^(X M ) = n 5 CV T 

2 Actually, the operator is only antisymmetric if the fermion is in a real representation, 
such as an adjoint, instead of a pseudoreal representation. Thus our method is consistent 
with Witten's result Jl6| that a theory of a single Weyl pseudoreal fermion is sick. 
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which implies the conventional (Euclidian) Majorana constraint on the 4-d 
fermion fields: 



Vv, W = Cffi( Xfl ) . (12) 

Using the same technique as in the Dirac case to remove bulk modes, we 
arrive at a formula for the pfaffian of a massless Majorana fermion: 

Pf \iK 5 C( p 4 + 75 <9 5 + Me(fl))] Pf [iC gj 

Pf\iK 5 C(p 4 + l5 d 5 + M)} Pf[iC(pt + M)) 1 ] 

This formula is easily extended to the lattice by replacing the Dirac action by 
the Wilson action in all five dimensions [|]]. As mentioned before, this leads 
to an answer identical to that derived by Neuberger [I2|, although derived 
in a somewhat different way. 

By using Neuberger's closed expression for the domain wall determinant, 
it is possible to show that the lattice version of the above pfaffian is positive 
definite, and hence can be computed unambiguously as the square root of the 
Dirac determinant []. Thus the domain wall approach has an added advantage 
over the Wilson fermion strategy, which suffers from a pfaffian which is not 
positive definite [I7|. it is therefore feasible with present technology to begin 
exploring this interesting theory. 



3 N = 1 SUSY Yang-Mills theory in d = 3 
dimensions 

N = 1 SYM in d = 3 is an interesting theory as well, especially in light 
of Witten's recent discussion of dynamical SUSY breaking ||18|| . Once again 
the spectrum consists of a gauge symmetry and a Majorana fermion, the 
gaugino. There are two independent relevant operators that break SUSY: 
the gaugino mass and the (quantized) Chern-Simons term, with one linear 
combination of the two being supersymmetric. In what follows we will assume 
that form some gauge groups it is possible to formulate the lattice theory such 
that the coefficient of the Chern-Simons term in the effective 3-d continuum 
theory vanishes (work in progress here!). In that case, the only relevant 
SUSY breaking operator is once again the gaugino mass. If we can realize 

3 This observation was made to DK by Y. Kikukawa. 
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chiral symmetry and gauge symmetry with a Majorana fermion, SUSY will 
once again arise in the continuum as an accidental symmetry, modulo the 
unresolved issue of the Chern-Simons term. 

We saw above that without constraints, a 5-d domain wall theory led to 
a massless Dirac fermion in 4-d; to end up with a Majorana fermion we had 
to impose a generalization of the Majorana constraint, which effectively took 
the square root of the 5-d domain wall determinant. However, following the 
same procedure in one fewer dimensions, a 4-d domain wall system with a 
Dirac fermion gives rise to two massless Dirac fermions in 3-d, four times as 
many degrees of freedom as we wish! In particular, 



where on the left hand side, the index % runs from 1 to 3 and the 7 matrices are 
4x4; on the right hand side, p 3 is the 3-d Dirac operator (2x2 dimensional in 
spinor space). Therefore it is clear we need to impose two binary constraints 
on the system. 

First of all, instead of using 4-d Dirac domain wall fermions, we can 
impose the 4-d Euclidian Majorana constraint, ip = C^ 7 \ where C 4 is a 4-d 
charge conjugation matrix. This naturally gives rise to a 3-d theory with two 
Majorana fermions localized at the two kinks. To reduce the spectrum to a 
single Majorana fermion in 3-d we use the trick of the previous section and 
constrain the field further to be Majorana under the 3-d charge conjugation 
matrix C3, and a simultaneous reflection TZ 4 in the compact fourth dimension. 
Thus the simultaneous constraints are: 

1. ^(xi,x 4 ) = C A ^ T (xi,x A ), 

2. ■%{x h x i )=K 4 C z $! T {x h x A ) 

To be explicit, one can choose the 7 matrix basis: 

7i = Ol <g) (J* , 74 = 3 <g) 1 , C 3 = 1 <g) 02 , C 4 = 0i <g> 02 . (15) 

It isn't obvious how to simultaneously impose these two constraints until 
one uses constraint (1) to replace constraint (2) by 



det [z(A7i + 74fl4 + Me(fl))] 
det [i{Dfii + 7 4 (9 4 + M)\ 



[detifl 3 ] 2 



(14) 



[deti(^ 3 + M)] 2 



2'. y(x h x 4 ) =TZ 4 C 3 C 4 1 ^(x l ,x 4 ) 
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This last constraint, relating \1/ to its reflection, tells us that we are living on 
an orbifold — only half the world we were considering represents independent 
degrees of freedom. So what we do is impose constraint (1) and compute the 
path integral over half our original space, namely for 6 e (0, it] with suitable 
boundary conditions at the fixed points of IZ^. 



1 - CsQ" 1 







(16) 



£4=0, ttR 



Then one finds the desired result, 
Pf [iC 4 (A7i + 7404 + Me(0))] 



Pf [iC 3 p 3 



Pf [zC 4 (A7i + 74<9 4 + M)\ Pf [zC 3 ( p 3 + M)\ 



6e(0,n] . (17) 



We have not yet completed our analysis of the reality/positivity of the 4-d 
pfaffians on the lattice, and the related issue of the Chern-Simons term in 
the effective 3-d theory. 



4 N = 2 SUSY Yang-Mills theory in d = 4 
dimensions 

N = 2 SYM in d = 4 would be fascinating to simulate on the lattice, since 
in the continuum it exhibits a vast array of interesting phenomena One 



might think that it impossible to do without fine tuning, however, because 
of the scalar fields in the N = 2 gauge multiplet. However, a promising idea 
is to formulate the theory first as an iV = 1 SUSY theory in d = 6 (starting 
from a domain wall theory in d — 7) The light spectrum of the d = 6 
theory, with UV cutoff A 6 would consist of gauge fields and a Weyl fermion. 
Then at a scale A 4 Ae, one compactifies to d = 4: the extra two gauge 
boson polarizations become the complex scalar of the d — 4, N = 2 gauge 
multiplet, while the Weyl fermion in d = 6 becomes the required two Weyl 
fermions in d = 4. Furthermore, all gauge, 4 and Yukawa couplings in the 
d = 4 effective theory are derived from the d = 6 gauge coupling g 6 . 

This approach is made respectable by the fact that in the continuum, the 
N — 1 SUSY algebra in d = 6 reduces under compactification to the N = 2 
SUSY algebra in d = 4 @. 



Of course, the idea is still to have the target N = 2 theory arise as an 
accidental symmetry in the effective theory. What one must try to do then 
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is to take A 4 sufficiently smaller than A 6 so that by the time one has scaled 
down to A 4 and passed over to the d = 4 effective theory, the theory is "su- 
persymmetric enough" to ensure that the noxious scalar masses radiatively 
generated in the effective d — 4 theory are "small enough" . 

How small is "small enough" ? To study the N = 2 theory in the strongly 
coupled region, where it is interesting, we need the scalar mass m s to satisfy 
m s <C A S qcd where Asqcd is the scale where the N = 2 gauge interactions 
get strong. 

Unfortunately this is impossible to achieve. The N — 1 supersymmetry 
in the d = 6 theory is only a symmetry of the operators of leading dimension; 
SUSY is violated by higher dimension operators, suppressed by powers of A 6 . 
Thus the SUSY violating radiatively generated scalar masses in the d = 4 
effective theory will be suppressed by powers of A^/ A§. We can suppress these 
terms as much as we want, by taking this ratio to be very small! However, 
the mass scale Asqcd is always smaller as it is exponentially small in A 4 /A 6 . 

To understand this, define the dimensionless gauge coupling g % = g 6 A 6 in 
the d = 6 theory. Since we begin with a weakly coupled domain wall fermion 
in d — 7 g§ ^ 1. The coupling of the d = 4 theory renormalized at the 
compactification scale A 4 is then given by g 4 = g& A 4 = g^A^/ Aq. Therefore 

Asqcd ~ A 4 e- 8 ^ ~ A 4e - 8 - 2 /^( A «/^) 2 « A 4 e~^ A ^ . (18) 

We see that while we obtain scalar masses suppressed by powers of A 4 /A 6 , 
the strong interaction scale Asqcd is exponentially suppressed in the same 
ratio. It follows that one cannot study the N = 2 theory in the interesting 
strongly interacting regime starting from a weakly coupled domain wall in 
d = 7, without fine tuning. 

The above argument does not rule out studying N = 2 SYM in d = 3 
by compactifying a d = 4 theory with approximate N — 1 supersymmetry, 
since the gauge coupling in d — 3 does not run logarithmically. However, 
this d — 3 theory has no ground state in the continuum, and so it does not 
seem interesting to simulate. 

5 Conclusions 

Domain wall fermions offer a compelling advantage over Wilson fermions in 
simulating N — 1 supersymmetric Yang-Mills theories on the lattice in d — 4 
and d = 3. In each case, supersymmetry arises as an accidental symmetry, 



9 



without fine-tuning. Both of these theories should be interesting to study in 
the near future. 

As for SUSY theories with scalars: it is hard to imagine how one can evade 
fine-tuning — after all, if one did have such a method, it would provide an 
alternative to SUSY as a solution to the hierarchy problem! 

It would be interesting to study perfect supersymmetric actions to try to 
extract the analogue of a Ginsparg- Wilson relation for supersymmetry, for 
then one might identify a clever approach to SUSY theories with scalars in 
the spectrum, one that minimizes the fine-tuning problems. 
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